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ABSTRACT 

Context. The gas mass of protoplanetary disks, and the gas-to-dust ratio, are two key elements driving the evolution of these disks 
and the formation of planetary system. 

Aims. We explore here to what extent CO (or its isotopologues) can be used as a tracer of gas mass. 

Methods. We use a detailed gas-grain chemical model and study the evolution of the disk composition, starting from a dense pre- 
stellar core composition. We explore a range of disk temperature profiles, cosmic rays ionization rates, and disk ages for a disk model 
representative of T Tauri stars. 

Results. At the high densities that prevail in disks, we find that, due to fast reactions on grain surfaces, CO can be converted to less 
volatile forms (principally s-COi, and to a lesser extent s-CH 4 ) instead of being evaporated over a wide range of temperature. The 
canonical gas-phase abundance of 1CT 4 is only reached above about 30 - 35 K. The dominant Carbon bearing entity depends on the 
temperature structure and age of the disk. The chemical evolution of CO is also sensitive to the cosmic rays ionization rate. Larger 
gas phase CO abundances are found in younger disks. Initial conditions, such as parent cloud age and density, have a limited impact. 
Conclusions. This study reveals that CO gas-phase abundance is heavily dependent on grain surface processes, which remain very 
incompletely understood so far. The strong dependence on dust temperature profile makes CO a poor tracer of the gas-phase content 
of disks. 
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1. Introduction 


Measuring the mass of protoplanetary disks remains a challeng¬ 
ing task. FF being either unobservable in the physical conditions 
prevailing throughout most of the disks, or when detected, only 
sampling peculiar conditions (Bary et al. 2008[ >. we have to rely 
on indirect tracers to estimate the disk masses. 


A common practice is based on the dust emission at mm 
wavelengths, where the dust opacity is moderate and thus the 
emission sensitive to most of the dust mass. However, the dust 
emissivity is not well known: it depends on dust characteris¬ 
tics, in particular the grain size distribution, and differential grain 
growth has been observed in disks (Guilloteau et al .[20i l||Perez[ 
|et al.|[20l2 !, invalidating the use of a unique emissivity coeffi¬ 
cient throughout the disk. Despite this limitation, under the sim¬ 
ple assumption of an average emissivity of order 2 cnrg -1 , it has 
been shown that the dust masses of protoplanetary disks scale 
with their host stellar mass (|Andrews et al.|20T3j). 


In any case, dust only cannot inform us about the gas con¬ 
tent. It is often assumed that the gas-to-dust ratio is of order 100 
to derive disk masses, but many processes (such as the vertical 
dust settling and the radial drift of large grains) may affect the 
gas-to-dust ratio. HD is perhaps the best proxy for Hi, and has 
been detected so far in TW Hya ( jBergin et al.|2013] >. However, 
the interpretation is strongly dependent on the thermal profile of 
the disk. Furthermore, there is no longer any sensitive enough 
telescope to detect it in more distant protoplanetary disks. 

It has been realized in the early studies of disk emission that 
the 4 CO lines were much fainter than expected from mm con¬ 
tinuum emission from dust (Beckwith et al. ||1990 ; Dutrey et al. 


1996). Part of this discrepancy may be attributable to different 


disk sizes in dust and continuum ( |Pietu et aL]2 ()06 2007; Hughes 
et al.|2009 ), but depletion of CO on dust grains was invoked as 


the main cause. 

Recently, Williams & Best ( |2014| suggested that, because of 
its relatively simple chemistry, CO could be a reliable tracer of 
dust mass. They argued that gas-phase CO abundance depends 
mostly on gas temperature, and that with a reasonable knowl¬ 
edge of the disk thermal structure, the fraction of carbon locked- 
up in gas-phase CO can be evaluated (10 4 where T > 20 K). 
Apart from the upper layers, which contain little mass, the ther¬ 
mal stmcture in disks is dominated by dust (Chiang & Goldreich| 


1997), and can be reasonably well constrained from SED mea¬ 
surements. Exploring a large grid of disk parameters, Williams 
|& Best| ( |2014| ) concluded that disk gas masses were in general 
rather small, leading to an average gas-to-dust ratio well below 
the canonical value of 100 found in the interstellar medium. 

However, this finding contradicts the high mass suggested 
by the analysis of the HD observations of TW Hya ( |Bergin[ 
et al. |201 3| >. To resolve the discrepancy between these two ap¬ 
proaches, a very low CO to H 2 ratio is required. Very low ap¬ 
parent CO to dust ratios were also found by Cha pillon et al.| 
(2008 1 for CQ Tau and MWC 758, two Herbig Ae stars, as well 
as for BP Tau ( |Dutrey et al.||2003] >, despite large temperatures 
(above 30 K) deduced from CO. Under the simple assumption 
that CO is mostly in the gas-phase because of the high temper¬ 
atures, ^F^IIIon£rar](j20T0| showed that the lack of detection 
of Cl in CQ Tau excludes photodissociation as the sole cause 
of low CO column densities, and implies a very low gas-to-dust 
ratio. 
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The CO rotational transitions being in general quite optically 
thick in disks, using the less abundant isotopologues is required 
to sample the gas content. This introduces a significant complex¬ 
ity because of the different shielding of the isotopologues against 
photodissociation by the UV radiation field. Another added com¬ 
plexity is the 13 CO fractionation which can occur in luke-warm 
regions (20-30 K). These processes have been studied in detail 
i Miotello et al.[p014|, who showed that these effect have sub¬ 


stantial impact on the derivation of the CO gas content, and could 
reduce the discrepancy between the two approaches. 

Favre et al. (2013) concluded that the emission of C 18 0 in 


TW Hya could only be explained assuming CO has been con¬ 
verted to other hydrocarbons, rather than released from grains. 
As an other alternative, jChapillon et al. ( 2008[ ) suggested that 
because large grains can remain cold, CO might remain trapped 
on such grains even in relatively warm disks, like those around 
the Herbig Ae stars CQ Tau or MWC 758. 

Given that other molecules suffer from an even more com¬ 
plex chemistry than CO and its isotopologues, understanding to 
what extent CO can be a tracer of disk masses is important. 

In this paper, we re-explore this issue using a detailed gas- 
grain chemical network. We concentrate here on some aspects 
of chemistry which can affect the (gas-phase) CO abundance in 
protoplanetary disks. We first study the impact of temperature 
profile and the role of grain surface chemistry for disks of differ¬ 
ent ages. Secondly we evaluate the importance of initial condi¬ 
tions on the time dependent disk chemistry. 

Our purpose is not to build complete disk models, but to re¬ 
veal the importance of these effects on the chemical evolution of 
protoplanetary disks. The paper is organized as follows. In sec¬ 
tion 2, we desccribe our chemical model and the physical disk 
structure used for this work. In Section 3, we present the effect of 
the disk vertical temperature profile on the chemistry of carbon¬ 
bearing species, as well as the impact of the cosmic rays ioniza¬ 
tion rate and the age of the disk. We also study the sensitivity 
of our model to the initial C/O ratio, the gas density and the age 
of the parent molecular cloud. Section 4 contains a discussion of 
our work. Finally, in Section 5 we present our conclusions about 
this work. 


2. Modeling 

2.1. Chemistry 


2.1.1. Model Description 


We used the Nautilus chemical model ( |Semenov et al.| |2010; 
Reboussin et al. 20141 which computes the abundance of species 


as a function of time in the gas-phase and at the surface of the 
grains. The model includes pure gas-phase chemistry with bi- 
molecular reactions, ionizations and dissociations by direct im¬ 
pact of cosmic ray particle and UV photons (interstellar FUV 
photons and secondary UV photons induced by cosmic ray/FD 
interactions). We assume the standard cosmic ray ionization rate 
^cr= 1-3x10 l7 s '. The dust grains are represented by spherical 
particles with a radius of 0.1 /mi, a density of 3 g.cm~ 3 , and are 
made of amorphous olivine. The gas and dust temperatures are 
assumed to be the same and we use a gas-to-dust mass ratio of 
100. The gas-grain interactions include adsorption of gas-phase 
species, those species can desorb back into the gas-phase by ther¬ 
mal desorption and non-thermal desorption (via the energy re¬ 
leased by exothermic surface reactions and cosmic ray induced 
desorption). The impact of photodesorption is discussed in Sect. 
3.2. The model takes into account grain-surface reactions with 


diffusion reactions at the surface of the dust particle via thermal 
hopping (we do not take into account diffusion by quantum tun¬ 
neling in this study). The energy barrier between two adjacent 
sites is taken as half of the binding energy (Garrod & Herbst 


2006). Photodissociation on the surface by UV photons and cos¬ 


mic ray induced UV field are also considered. 


2.1.2. Cloud Model Parameters 


To obtain the initial chemical composition of the disk, we 
first compute the chemical composition of the parent molec¬ 
ular cloud. We run Nautilus during 10 6 yr for typical dense 
cloud conditions as our nominal cloud: a gas density of 
n H = n(H) + 2n(FD) = 2xl0 4 cm' 3 , a temperature of 10 K (with 
7'„ as = r dust ), a visual extinction of 10, a cosmic ray ionization 
rate of 1.3x10 17 s 1 and the initial abundances as listed in Table 
T] For all elements but He, N, C and O, we have used the “low 
metal” elemental abundances defined by Graedel et al. ( |1982| . 
For oxygen, two values have been considered: a “low depletion” 
case which corresponds to C/O = 0.5 in our standard model, and 
a “high depletion” case which corresponds to C/O =1.2 consid¬ 
ered in Sect. 13.5.21 The helium abundance is assumed to be 0.09 
( Wakelam & Herbst|2008| ) while the nitrogen and carbon abun¬ 
dances are extrapolated from Jenkins ( 200 9) up to a density of 
2xl0 4 cm' 3 (see discussion in Hincelin et al. 2011). 


Table 1. Initial abundances used in our model. 


Element 

Abundance relative to H 

H 2 

0.5 


He 

9x10' 2 R 


N 

6.2x1 0^ 

b 

O 

3.3xl0~ 4 

(C/O=0.5) 


1.4xl0~ 4 

3 (C/0=1.2) 

C + 

1.7xl0- 4 @ 

s + 

8x10~t] 


Si + 

8xl0' 9 ® 


Fe + 

3x10' 9 E 


Na + 

2x1 O' 9 


Mg + 

7x10' 9 >3 


P + 

2x10' 10 I3 

Cl + 

1x10' 9 >3 



See discussion in 


Jenkins 


(2009). 


See discussion in 


Wakelam & Herbst 

2008 


Hincelin et al. 

2011 

). 


d Low metal abundances (Graedel et al. 


1982). 


2.1.3. Chemical Network 

The chemical network used for this work contains 8624 reac¬ 
tions: 6844 are pure gas-phase reactions and 1780 are grain- 
surface and gas-grain interactions. The model follows the chem¬ 
istry of 703 species (atoms, radicals, ions and molecules): 504 
are gas-phase species and 199 are species on grains. The surface 
network is based on |Garrod et al.| ( |2007| ) whereas the gas-phase 
network is based on kida.uva.2011 with updates from Wakelam 
[et al.| ( |20131 ) and |Loison et al.| ( ]2014a|b) i. We have added 19 
gas-phase species and 377 pure gas-phase reactions (mainly car- 
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bon and nitrogen chains). The network is available on the KIDA 
(Kinetic Database for Astrochemistry) website [j] 


2.2. Disk Structure 

For the disk physical parameters, we used the disk model de¬ 
scribed in Hersant et al. ( |20091 >. 

The vertical density structure is obtained by integrating the 
equation of hydrostatic equilibrium: 


( 1 ) 


where G is the gravitationnal constant, M, = 0.5M 0 is the mass 
of the central star, // = 2.4 is the reduced mass of total hydrogen, 
bih the mass of the proton and k B the Boltzmann constant. 

For the isothermal case, the vertical density is given by: 


d In p 

(GM„z\ 

A"«H 'j 

, 51n7 1 

dz 

l r 3 ) 

Ub t) 

dz 


pfe)=p„exp-(iy 


( 2 ) 


where 77 = feT^r 3 /GM t /umu) 1 ^ 2 is the pressure scale height 
derived from the midplane temperature r m id. 

We assumed the disk surface density to vary as r~ ! 5 as de¬ 
scribed in |Hersant et al.| ( [2009} >: 


I(r) = I 


too 


(100au) 


-1.5 


gem 


(3) 


with Sioo = 0.8 g cm~ 2 . The hydrogen surface density is given by 
£ioo = lumnNioo(H) in which n is the mean weight per H nuclei, 
which we approximate as /r/2. Thus Num(H) s 4 x 10 23 cm~ 2 . 
With an outer radius of 700 au, the disk mass is = 0.03 M a . 

Contrary to Hersant et al. (|2009|l, we use th e disk temper¬ 
ature structure as described in [Williams & Best| ( j2014j i. The ra¬ 
dial midplane temperature 7' m j c i and the atmospheric temperature 
r atm (temperature at z > 477) are given as a power law of the ra¬ 
dius: 

T m id(r) = Tmidj (t^—) (4) 

V1 au/ 


T’atmM = T atm , -- 


(i=)' 


(5) 


with q = 0.55, which leads to a flaring exponent of 1.23. 

A sine function is used for the connection between the mid¬ 
plane and the atmosphere gas temperatures: 


T(r,z) 


1 T m id + (T’atm - 7’ m id) 

sin (f) 

[ T atm 



26 


ifz<z q (6) 

if z>z q 


with <5 = 2 (steepness of the profile) and z q = 4// (height at 
which the disk reaches the atmospheric value). The visual ex¬ 
tinction is related to the total column density of hydrogen by 
A v = N(H)/ 1.6 x 10 21 cm~ 2 , following Wagenblast & Hartquist 
( |1989| >, consistent with the adopted mean grain radius of 0.1 jjm. 
Changing the temperature profile also affects the density and vi¬ 
sual extinction distributions. 

For the UV radiative transfer in the disk, we consider only 
the vertical extinction. We assume that the radial extinction of 
stellar UV field is very efficient so that UV photons are scattered 
by small grains at high altitudes. The UV flux is decreasing as 
1/r 2 from the central star. We consider Go = 410 for the field 
for a distance of 100 au from the star: this is representative of 





http://kida.obs.u-bordeaux 1 .fr/models 


Av 

Fig. 1. Examples of temperature, density and visual extinction 
vertical profiles in the protoplanetary disk at 300 au used for 
the chemical modeling. Solid lines are for Model 1 and dashed 
lines for Model 2 (see Table [2ji. The vertical profiles for these 
two models are shown for three different values of the midplane 
temperature 7) ni( j at 300 au: ~ 4 K (blue lines), ~ 20 K (green 
lines) and ~ 40 K (red lines). Note that the disk heights depend 
on the temperature: from the coldest to the warmest model, the 
disk heights are respectively 29, 62 and 88 au at 300 au. 


low mass T Tauri stars like DM Tau. Half of this flux is assumed 
to be scattered downwards inside the disk. To model the CO and 
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We created a large grid of models by varying the values of 
7 mK ij and 7 atm j. The disk midplane temperature is therefore 
ranging between 4 and 52 K while the temperature in the up¬ 
per layer is ranging between 17 and 95 K (at 300 au). We also 
considered different models (see Table [2] for a description of our 
models). Model 1 is considered as our nominal model. In Model 
2 we study the impact of the steepness of the temperature profile 
taken as stronger in this model by increasing the value of T’atm.l- 
A third model (Model 3) was considered in which the dust tem¬ 
perature has been considered as different from the gas one. In 
that case, the gas temperature is assumed to follow the strong 
gradient temperature of Model 2 while the dust temperature fol¬ 
lows the one of Model 1. Fig|T] shows the disk temperature, vi¬ 
sual extinction and density vertical profiles at 300 au for some of 
those models. Depending on the temperature profile considered, 
disks have different height, the warmer ones are thicker. From 
the coldest to the warmest one shown in Fig|T] the disk heights 
at 300 au are respectively, 29, 62 and 88 au. 


3. Results 

In this section, we present the results obtained with our models. 
We first study the impact of the temperature profile on the disk 
chemistry and more particularly on the carbon-bearing species 
abundances. Secondly, we show the sensitivity to the initial 
cloud composition. The initial chemical conditions used for the 
disk modeling are listed in Table [4] Hereafter, for any molecule 
X, we shall use s-X to refer to molecules trapped on grains, and 
X for molecule in the gas-phase. 


3.1. Impact of the vertical temperature profile 

The (gas-phase) CO vertical abundance at r = 300 au is repre¬ 
sented in Fig|2]for the 6 profiles shown in FigjT|and for a disk 
age of 1 Myr. In the upper layer of the disk, the CO abundance is 
mostly affected by the photodissociation process. Because of the 
different scale heights, the visual extinction at Z/H -4 is vary¬ 
ing between models (see Fig. |T]». In the upper layer, CO is more 
efficiently photodissociated by UV photons for the warmer mod¬ 
els (green and red ones) because the visual extinction is smaller. 
The largest CO abundances are obtained in the warm molecular 
layer, whose height and thickness differs according to the model. 
The coldest model with a stronger temperature gradient (blue 
dashed line in Fig|2]> shows a larger CO abundance at smaller 
scale height compared to the coldest model with a weak tempera¬ 
ture gradient (blue solid line in Fig(2]). This is due to the fact that 


Gas-Phase CO abundances 

Fig. 2. Vertical gas-phase CO abundances (relative to total hy¬ 
drogen) at 300 au for a disk age of 1 Myr. The color coding 
corresponds to the vertical profiles show in Fig|T| Solid lines are 
for Model 1 and dashed lines for Model 2 (see Table [2]). 


the desorption temperature of CO is reached at smaller height. 
In the disk midplane, the larger CO abundance is obtained with 
the higher T ml( j model (in red in Fig JTJ because of a more effi¬ 
cient desorption of s-CO from the ices. For this model, the CO 
depletion is small, as expected at such high temperatures. The 
CO abundance obtained with the smaller T mllJ model (in blue in 
Fig jT]> is higher than the intermediate model (green one). In fact, 
around 15 K, s-CO is destroyed on the dust grains to form more 
complex species by secondary UV photons induced by cosmic- 
rays, such as S-CO2, s-HiCO and S-CH3OH, whereas when the 
temperature is too low (below 9 K), the diffusion of the species, 
even atomic hydrogen, at the surface of the grains is less effi¬ 
cient. 

The CO abundance profiles obtained with the models are 
then sensitive to the vertical temperature profile that is assumed: 
the midplane temperature but also the steepness of the gradient. 
For midplane temperatures larger than 20 K, the steepness of 
the gradient no longer substantially affects the CO abundance 
because the evaporation temperature of s-CO is reached in the 
molecular layer. 

For the models displayed in Fig. [2] only the warmest ones 
(red lines) reach the value CO/H 2 = 10 4 . For the other models, 
which are colder, this ratio is much smaller for two main 
reasons. Firstly, the warm molecular layer in which larger CO 
abundances are obtained is less extensive than in the case of 
warmer models and secondly, the CO depletion due to freeze-out 
and grain surface chemistry is very efficient in the disk midplane. 

Figure [3] represents the disk column densities and vertical 
abundances (computed from the molecular column densities and 
the total H column density) of some carbon-bearing species in 
the gas-phase and at the surface of the grains as a function of the 
disk midplane temperature and for the three models as described 
in Table [2] The results are obtained for a disk age of 1 Myr. 

In the case of temperature profiles with a midplane tempera¬ 
ture between 5 and ~ 15 K in Model 1, s-CO is the most abun¬ 
dant C-bearing species. At higher temperature, its abundance de¬ 
creases because it is converted into other C-bearing molecules 
such as S-CH 4 . The direct photodissociation of CO in the gas 


4 















Reboussin et al.: Chemistry in Protoplanetary Disks 


Gas-Phase Species 



Surface Species 



Fig. 3. Average abundances and integrated column densities 
(cm -2 ) of carbon-bearing species in the gas-phase and on the 
grain surfaces for Model 1, 2 and 3 as a function of the disk 
midplane temperature (T m jd). s- refers for surface species. Those 
results are obtained at 300 au for our nominal cloud and for a 
disk age of 1 Myr. 

and at the surface of the grains increases the abundance of neu¬ 
tral carbon, and subsequent hydrogenation of Carbon on the sur¬ 
faces contributes to form S-CH 4 . Above 15 K, S-CO 2 becomes 
the most abundant molecule. Its formation is due to the follow¬ 
ing reactions, for which the efficiencies are increased by 10 to ~ 
20 orders of magnitude at such temperatures compared to 15 K. 

s-O + s-CO —> S-CO 2 (7) 

s-O + s-HCO -> S-CO 2 + s-H ( 8 ) 

s-OH + s-CO -> S-CO 2 + s-H (9) 

Reaction[7]is the main formation path of S-CO 2 in our model al¬ 
though it requires to overcome a large activation energy barrier 
(1000 K). Reaction [ 8 ] being barrierless contributes also to form 
S-CO 2 efficiently. Reaction [9] requires to overcome a small acti¬ 
vation energy barrier (80 K). The warm dust temperature allows 
to increase the mobility of species at the surface of the grains 
and the diffusion becomes much more efficient. However, the 
dust temperature is not warm enough to desorb the S-CO 2 back 


Table 3. Molecular binding (desorption) energies. 


Species 

E D ( K) 

C 

800 

CO 

1150 

ch 4 

1300 

h 2 co 

2050 

co 2 

2575 

CH 3 CCH 

3837 

C 2 H 6 

4387 

CH 30 H 

5534 


into the gas-phase. At even higher temperatures (above 30 K), 
CO in gas-phase is the reservoir of C because of efficient des¬ 
orption and the lack of sticking on the grains. 

Model 3 (in which T gas > 7j ust ) produces species abun¬ 
dances very similar to Model 1 whatever T m ;d, meaning that it 
is the grain temperature that is the important parameter here. For 
temperature between 25 and 50 K (depending on the molecule), 
the steepness of the temperature gradient does not influence the 
computed abundances in the gas-phase or at the surface of the 
grains. For T m ;d smaller than 25 K, the gas-phase abundances 
are larger for Model 2 mainly because s-CO desorbs from the 
grains in Model 2 in the molecular layer, producing the other 
C-bearing species, while it stays on the grains in Model 1. 

S-CH 3 OH and S-H 2 CO are efficiently formed at low T m jd: 
between 10 and 20 K, their abundances are respectively, ~ 10 1 
and ~ 10 6 . Their main formation paths are: 


s-OH + s-CH, -> S-CH 3 OH 

( 10 ) 

s-H + s-CH 2 OH -> S-CH 3 OH 

( 11 ) 

s-H + s-HCO -> s-H 2 CO 

( 12 ) 

s -0 + S-CH 2 —* S-H 2 CO 

(13) 


However, their abundances are strongly decreased at higher T m ;d 
since s-H, s-HCO, s-O and S-CH 2 are easily desorbed at such 
temperatures. The largest CH 3 OH and H 2 CO gas-phase column 
densities are found to be around 10 12 cm -2 for T m jd = 20 K. 

In summary, the carbon reservoir in protoplanetary disks 
strongly depends on the disk vertical temperature profile which 
is controlled by the midplane temperature; the reservoirs are: 

- s-CO in the ices below 15 K. 

- S-CO 2 in the ices between 15 and 30 K. 

- CO in the gas-phase above 30 K. 

The results described above were obtained at 300 au. We 
show now the results of Model 1 (weak temperature gradient) 
at 100 au. In Fig J4] are represented the disk abundances and col¬ 
umn densities of carbon-bearing species in the gas-phase and on 
the grain surfaces at 100 au for Model 1 as a function of the disk 
midplane temperature (to be compared to Fig|3j. At 100 au, for a 
given vertical temperature profile, the density, the visual extinc¬ 
tion, and the UV fluxes are larger than at 300 au. Considering 
the larger densities, the reservoirs of carbon at low temperature 
(T m id < 17 K) are S-CH 4 and s-CO, and S-CO 2 between 17 and 
34 K. At larger temperatures, gas-phase CO becomes dominant. 
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Disk Midplane Temperature (K) 


Fig. 4. Average abundances and integrated column densities 
(cm -2 ) of carbon-bearing species in the gas-phase and on the 
grain surfaces at 100 au as a function of the disk midplane tem¬ 
perature (Tmid). Results are for Model 1 and for a disk age of 1 
Myr. 


3.2. Photodesorption 


Photodesorption was not included in our previous models. High 
desorption yields were found by Oberg et al. (2007), but it has 
been shown that the desorption efficiency depends on the shape 
of the UV spectrum (Fayolle et al. |201 1 ). In particular, the pho¬ 
todesorption yield of CO ice induced by Ly-cr photons is small 
( jChen et al.|20l4| . 

We tested the impact of photodesorption on our results at 100 
and 300 au for a disk age of 1 Myr. We used a large yield of 10 5 
molecules photons -1 for all species except for CH3OH, HiO, 
CO7, N?, and CO for w hich we used the formalism described in 
Oberg et al. (2009a b|c I. 

The carbon reservoir remains the same for both radii. The 
gas-phase CO abundance is only affected at 300 au for midplane 
temperatures between 4 and 15 K, with a CO abundance five 
times larger than in our nominal model. However, s-CO remains 
the main carbon-bearing species. For both radii, only more com¬ 
plex molecules, such as CjHft and CH 3 CCH, are efficiently pho- 
todesorbed (their abundances are increased by more than a factor 
10 ). 


3.3. Impact of the cosmic ray ionization rate 


Cosmic rays are important for the ionization in the disk. 
However, as recently suggested by Cleeves et al. ( |2013) l, cosmic 
rays penetration may be highly reduced in the presence of stellar 
winds and/or magnetic fields. In this section we study the impact 
of a lower (cosmic ray) ionization rate on the disk chemistry. 

Figure [5] shows the abundances of the main Carbon-bearing 
species for two values of the cosmic rays ionization rate: 
1.3xl0 -17 and 1.3xl0 -18 s -1 (respectively in solid and dashed 
lines). The impact of ionization rate on CO abundance is limited, 
being restricted to the temperature range 20 - 35 K. Considering 
the lower value for ^cr. s-CO abundance increases because its 
dissociation on the grain surface by UV photons induced by cos¬ 
mic rays is less efficient. As a consequence, there is less carbon 
available to form S-CH 4 . In the gas-phase, the dissociation by 
secondary UV photons induced by cosmic ray/H 2 interactions is 
also less efficient so that CH 4 abundance increases at 300 au. At 
100 au, the UV flux being larger, the direct photodissociation is 
the dominant process. CH 4 presents therefore similar sensitivity 
than S-CH 4 and decreases as well. 


We also ran our Model 1 for £ C r = 1-3x10 


19 


but the 


reservoirs remain the same as those shown in Fig|5] (dashed 
lines) for both disk radii. 


3.4. Impact of the age of the disk 

The conversion of s-CO into more complex molecules on 
grains is a time dependent process, which requires a minimum 
timescale to occur. Up to now, the results described were ob¬ 
tained for disks of 1 Myr. In this section we study the chemistry 
of carbon-bearing species for disks 5 times older. 

Figure [ 6 ] shows the abundances of CO, s-CO, S-CO 2 , and s- 
CH 4 , obtained for Model 1 and for Tmid = 20 K, as a function 
of time. The s-CO abundance strongly decreases at 10 s yr due to 
its conversion into other species such as S-CO 2 and S-CH 4 . We 
can see that the formation of S-CH 4 requires a longer timescale 
to occur than the formation of S-CO 2 . S-CH 4 reaches its highest 
abundance at 10 7 yr while S-CO 2 is the most abundant species 
for a disk at age of 10 6 yr. The reservoirs of carbon are therefore 
changed while considering an older disk. 

Figure [7] represents the abundances and column densities of 
the main carbon-bearing species as a function of the disk mid¬ 
plane temperature for disks of 1 and 5 Myr. For the older disk, 
S-CH 4 is now the reservoir of carbon between 10 and 22 K in the 
midplane which is due to its long formation timescale. S-CO 2 is 
also more abundant than CO between 32 and 40 K at 5 Myr com¬ 
pared to 1 Myr. This apparent dependency with temperature is 
actually a density effect. Because of the hydrostatic scale height, 
densities in the disk are smaller for higher temperatures, and the 
conversion time is just inversely proportional to the density. The 
conversion of CO into S-CO 2 is then longer in warmer disk so 
that at 35 K CO is the main carbon-bearing species for a disk 
of 1 Myr. For instance, in our model, this conversion timescale 
varies from 6 x 10 5 yr for disks at 25 K to 2 x 10 6 yr for disks at 
35 K. 


3.5. Sensitivity to the initial conditions 

We now study the impact of the chemical composition of the par¬ 
ent molecular cloud on the disk chemistry. Several parameters 
have been studied such as gas density, C/O elemental ratio and 
age of the parent cloud. The abundances of the main C-bearing 
species obtained for different cloud parameters, and used as ini- 
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300 AU 



100 AU 



Fig. 5. Average abundances and integrated column densities 
(cm -2 ) of main reservoirs of carbon at 300 au and 100 au for 
two values of cosmic-rays ionization rate as a function of the 
disk midplane temperature (T mld ). Results are for Model 1 and 
for a disk age of 1 Myr. 



Fig. 6. Average abundances and integrated column densities 
(cm -2 ) of main reservoirs of carbon at 300 au as a function of 
time. The results are obtained for Model 1 and for T m [ d - 20 K 
at 300 au. 



Fig. 7. Average abundances and integrated column densities 
(cm -2 ) of main reservoirs of carbon at 300 au for two disk ages 
as a function of the disk midplane temperature (Tmid). Results 
are obtained for our Model 1. 


tial conditions for the disk chemistry, are listed in Table [4] The 
results presented in this section are for disks of 1 Myr. 

3.5.1. Density of the parent molecular cloud 

To test the importance of the density of the parent molecular 
cloud, we redid our simulations but multiplied the density of the 
parent cloud by a factor of ten. The disk carbon chemistry is not 
very sensitive to the gas density of the parent molecular cloud. 
The reservoirs of carbon as described in the previous section are 
not changed. Only CH 3 CCH abundance in the ices is decreased 
whatever T mld by a factor up to ~ 2. 


3.5.2. C/O elemental ratio 


For the carbon abundance, we use the value determined by 
|Jenkins ( |2009) > and we vary the C/O ratio by considering 
two different values for the oxygen elemental abundance: 1 ) 
3.3x10 4 (value of our nominal model), a low depletion case. 


and 2) 1.4x10 4 , a high depletion case (see section 2.1.2). 

The reservoirs of carbon are not changed while using the 
high depletion case for C/O elemental ratio (C/O = 1.2). The 
only main difference is a smaller abundance of CO containing 
species and a larger abundance of CH containing species under 
some conditions. At 7 K for instance, the CO, CH 4 and CO 2 
abundances in the ices are respectively 8.8x10 5 , 2.5x10 2 and 
3.8xl0- 7 . 


3.5.3. Age of the parent molecular cloud 

The age of molecular clouds is still a subject of debate. Ages 
are often considered between a few million years ( |Hartmarin| 
|et al.|200l]|Ballesteros-Paredes & Hartmann|2007|l and 10 mil- 
lion years (Mouschovias et al.J 2006j >. |Ballesteros-Paredes &| 
|Hartmann| ( |2007| ) suggested a short timescale of a few million 
years. They argued that starless molecular clouds being rarely 
observed, star formation should occur shortly after cloud for- 
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Gas-Phase Species 



.i.i.i.i.i.n 1 

0 10 20 30 40 50 60 


Disk Midplane Temperature (K) 




Disk Midplane Temperature (K) 


Fig. 8. Average abundances and integrated column densities 
(cm -2 ) of carbon-bearing species in the gas-phase and on the 
grain surfaces for three ages of the parent molecular cloud as a 
function of the disk midplane temperature. Results are for Model 
1 at 300 au, and for a disk age of 1 Myr. The results obtained for 
C and C + are not displayed since the differences are very small. 


mation. An argument favouring long cloud lifetime is the dif¬ 
ficulty to explain the sufficient formation of molecular hydro¬ 
gen in only a few Myr considering the slow rate at which FD 
forms in the ISM ( |Hollenbach et al.|1971[ l. However, |Glover &| 
Mac Low ( 2007[> suggested that turbulence significantlly reduces 


the time required for the conversion of H into H 2 (less than 1-2 
Myr). Pagani et al. (2011 1 suggested that the absence of DCO + in 
dark clouds is due to a high ortho-H 2 abundance which provides 
an upper limit to the age of 3-6 million years. More recently. 


Britnken et al. (2014) used the abundance ratio of ortho-H 2 D + 


and para-H 2 D + and determined an age of dense molecular clouds 
around one million year. 

To study the sensitivity of the model results to the age of 
the parent molecular cloud, we have considered two additional 
cases. We computed the initial cloud composition at 10 s vr and 
10 7 yr, in addition to our standard model at 10 6 yr. In Fig|8] we 
display the abundances and column densities of carbon-bearing 
species as a function of the disk midplane temperature computed 
with the different initial compositions. The abundances of the 
main C-bearing species obtained in the parent cloud at the 
different times (10 s , 10 6 and 10 7 yr) are listed in Table [ 4 ] The 
reservoirs of carbon in the disk are the same if a younger parent 
cloud is considered. However CO reaches its canonical value 
at larger 7" m id: 40 K instead of 32 K. In the case of an older 
cloud the reservoirs of carbon are changed. At temperatures 
smaller than 9 K for instance and in the case of a 10 7 yr cloud, 
CO in the ices is not the main C-bearing molecule anymore 
while CH 4 in the ices is. This can be explained by differences 
in the initial chemical composition of the cloud. The 10 s and 

10 6 yr compositions are similar and only the fraction between 
gas-phase and surface CO is reversed. At the large density of the 
disks, the gas-phase CO is quickly depleted on the grains. At 

10 7 yr in the cloud, the CO in the ices is hydrogenated to form 
H 2 CO and CH 3 OH in the ices but is also photodissociated by 
secondary UV photons induced by cosmic-rays and the atomic 
carbon is hydrogenated to S-CH 4 . So that the reservoir of carbon 
in the disk remains CH 4 in the ices. At larger temperature, 
S-CO 2 is not anymore the reservoir of carbon since CO 2 is 
formed on the surfaces from CO (already under abundant in the 
cloud) and the reservoir is C 2 H 6 in the ices. As a summary, the 
reservoirs of carbon in a disk formed from an older molecular 
cloud with an age of 10 7 yr are CH 4 in the ices for temperature 
profiles computed with a T lnld < 22 K and C 2 H 6 in the ices for 
larger T mid . 


The reservoir of carbon in disks is thus rather insensitive to 
initial conditions. The age of the cloud has an effect only if it is 
old enough, at least 10 Myr, a value on the high end of estimated 
cloud ages. In that case, s-CH 4 becomes the dominant C-bearing 
molecule instead of s-CO in the ices (for T ml( j < 22 K). For in¬ 
stance, for T m id = 15 K, the S-CH 4 /S-CO ratio is ~ 5. Such a ra¬ 
tio is much larger than those found in Solar System comets (see 
Bockelee-Morvan et al. 2004| ), which may indicate 


Table 1 of 
that our So 


ar System originates from a younger cloud. 


4. Discussion 

Transformation of CO into more complex species on grain sur¬ 
faces was already invoked by jAikawa et al.| ( |1997) >. However, 
in their study, the process is driven by neutral carbon and re¬ 
quired enhanced X-ray ionization to be started. Ionized helium 
atoms react with CO to release C, which further react on grains 
to get CO 2 or other hydrocarbons. We have shown here that with 
a more complete grain chemistry, this mechanism is no longer 
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necessary, and CO can be converted to other forms even without 
enhanced ionization. 

Our modeling reveals that, because of this processing on 
dust grains, the CO chemistry does not behave in a simple way. 
Figure [9] compares our average CO abundances with those pre¬ 
dicted using the simple prediction of |Williams & Best] ( |2014[ >, 
who assumed CO/H 2 = 0 where T < 20 K and in the pho¬ 
todissociation region, and CO/Hi = 10 4 elsewhere. Except for 
the coldest disks (where it predicts no CO), the simple prescrip¬ 
tion clearly overpredicts the average CO abundance for midplane 
temperatures below 30 to 35 K. At 20 K, the difference between 
our chemical model and the simple prescription is a factor of 
~ 40. This suggests that gas-phase CO abundance is critically 
dependent on the (dust) temperature distribution and grain sur¬ 
face chemistry. Figure Kljshows that an abundance of 10 4 is only 
reached for disks with a midplane dust temperature above 30 K 
at 1 Myr and above 40 K at 5 Myr. 


Gas-Phase CO 



Disk Midplane Temperature (K) 

Fig. 9. Average abundance of gas-phase CO as a function of the 
disk midplane temperature (T nmi ). The results obtained with our 
Model 1 for disks of 1 and 5 Myr and at 100 and 300 au are 
compared to the model used by Williams & Best;( ]2014) for CO 
chemistry. 


This complex CO chemistry may explain the low C lfi O to 
dust ratio observed in TW Hya by Favre et al. ( 2013|> . With 
an age estimate of 7- 8 Myrs ( |Barrado Y Navascues| |2006| 


Ducourant et al. 


2014 ) 2 ] TW Hya is old enough compared to 


the timescale required to convert CO to other forms on grain 
surfaces. To estimate quantitatively the fraction of CO in gas 
phase, a very accurate thermal model will be needed. |Hughes| 
et al. ( 2008} ) reproduced the mm continuum and CO images of 
the TW Hya disk by two different disk models, a “cold” one 
with T(r ) = 30K(r/100au) 05 , and a “warm” solution with 
T(r) = 40K(r/100au) 112 . We also note that Hughes et al. (2008 1 


assume 7’,| ust = 7 gils , and the temperatures are certainly biased 
towards somewhat higher values because of the expected verti¬ 
cal temperature gradient and of the difference between gas and 
dust temperatures. Nevertheless, for both cases, we would pre¬ 
dict abundances of order 10~ 5 , in reasonable agreement with the 


Vacca & Sandell 


( 2011 ). 


• a younger age of 3 Myr has been quoted by 1 
but seems inconsistent with the dynamical age derived by |Ducourant| 
etal. (20T4I) 


results of Favre et al. ( |2013) >. Thus, there may be no need to in¬ 
voke other mechanisms than the conversion of s-CO to s-CCB 
which is the dominant process in our study. 

When comparing for other objects, CO depletion due to sur¬ 
face chemistry may still have difficulties in explaining the appar¬ 
ent low CO to dust ratio derived by Chapillon et al. ( 2008| > for 
the very warm disk of CQ Tau, where temperatures may exceed 
50 K. However, Chapillon et al. ( 2008) only marginally resolved 
the source and their quoted temperatures are for gas only. Higher 
resolution studies, coupled to a more detailed thermal and chem¬ 
ical modeling for this source, would be needed to conclude. 

We did not consider the very strong UV flux of Herbig Ae 
stars, so that a direct comparison between our results and any 
specific object is premature. Similarly, we only considered small 
grains in this study. Larger grains will reduce the gas-grain in¬ 
teractions, thereby augmenting the conversion timescale for CO, 
but as pointed out by Chapillon et al. ( 2008) , these large grains 
could remain cooler. Also, turbulent mixing may play a role 
in bringing more frequently dust in the warm region, reducing 


the CO deficit. Turbulent mixing has been studied by Furuya 
& Aikawa (j2014|), and found to suppress the carbon depletion 


only within the inner 50 au. High CO content is thus expected 
in the inner regions (radius perhaps up to 30-50 au only). This 
is consistent with the CO/Hi ratio about 10 4 derived from IR 
measurements by France et akj(|2014) for RW Aur. 


5. Conclusions 

In this paper, we have discussed the chemistry of carbon-bearing 
species in protoplanetary disks using our gas-grain chemical 
model Nautilus. In particular, we studied the effect of the disk 
temperature structure on the chemistry of carbon reservoir in 
disks. We also studied the impact of different initial parent cloud 
compositions. We found that the reservoir of carbon in disks is 
not sensitive to initial conditions such as gas density or C/O ele¬ 
mental ratio of the parent molecular cloud. 

However, carbon chemistry is strongly sensitive to the disk 
vertical temperature profile. In many of our models, we have 
found that the CO gas-phase abundance is lower than the canoni¬ 
cal value of ~ 10 4 generally assumed. Indeed, with our models, 
such abundances are only reached for warm disks with a mid¬ 
plane temperature above 30 K, at 100 and 300 au for a disk age 
of 1 and 5 Myr. For smaller temperatures, we found that CO is 
converted into other species at the surface of the grains, such 
as S-CO 2 , for which higher freezeout temperatures prevent their 
desorption back into the gas-phase. The effect is slightly less effi¬ 
cient at low cosmic-rays ionization rate (down to 1.3x10 lx s 1 ) 
and large radii (300 au). 

Finally, our work demonstrates that CO is not a simple tracer 
of H 2 in protoplanetary disks, because its abundance is very sen¬ 
sitive to the large (dust) temperature gradients that exists there. 
These gradients will be very difficult to constrain with sufficient 
accuracy. Given the additional complexity of isotopologue selec¬ 
tive photodissociation and fractionation, evaluating disk masses 
from measurements of the CO isotopologue emission only will 
remain a very challenging task. 
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Table 4. Predicted abundances in the molecular cloud considering different model 
parameters. 


Abundances 0 

(with respect to n H ) 

species 

nominal clout 

C/0=1.2 

10 3 yrs 

10' yrs 

n H =2xl0 3 cm 3 

CO 

2.4(-5) 

1.3(-5) 

7.21-5) 

2.91-7) 

2.2(-6) 

ch 4 

1-5 (-7) 

4.0(-7) 

4.51-7) 

6.51-7) 

2.41-8) 

CO. 

3.3(-8) 

6.0(-9) 

4.71-8) 

9.2(-ll) 

7.91-10) 

c 

3.1(-8) 

7.0(-8) 

2.91-5) 

5.41-8) 

1.41-8) 

c + 

4.2(-9) 

8.0(-10) 

2.01-9) 

1.21-13) 

4.41-10) 

C 2 H 6 

6.2(-9) 

7-1(-9) 

7.11-11) 

1.21-8) 

8.31-9) 

h 2 co 

9.8(-9) 

1.8(-8) 

4.51-8) 

1.61-9) 

2.31-9) 

CH,OH 

7.1(-10) 

1.0(-9) 

7.21-11) 

6.21-11) 

1.11-10) 

s-CO 

7.2(-5) 

3.0(-5) 

2.51-5) 

1.91-7) 

9.91-5) 

s-CH 4 

8.5(-6) 

2.51-5) 

4.2(-6) 

9.21-5) 

5.11-6) 

s-CO. 

1.6(-6) 

5.61-8) 

8.11-7) 

3.91-9) 

2.31-6) 

s-C 2 H 6 

3.4(-6) 

3.91-6) 

3.91-8) 

6.4(-6) 

4.6(-6) 

s-H.CO 

2.8(-5) 

4.K-5) 

1.81-6) 

1.01-6) 

2.31-5) 

s-CH 3 OH 

9.7(-6) 

1.61-5) 

5.0(-7) 

2.51-6) 

5.81-6) 

S-CH3CCH 

2.4(-6) 

6.11-6) 

1.0(-6) 

4.61-6) 

1.01-6) 


a a(b) represents ax 10* 
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